A WOX homolog disrupted by a transposon led to the loss of spines and contributed to the domestication of lettuce.

The loss of spines is one of the most important domestication traits for lettuce (Lactuca sativa). However, the genetics and regulation of spine development in lettuce remain unclear. We examined the genetics of spines in lettuce using a segregating population derived from a cross between cultivated and wild lettuce (Lactuca serriola). A gene encoding WUSCHEL-related homeobox transcription factor, named as WOX-SPINE1 (WS1), was identified as the candidate gene controlling the spine development in lettuce, and its function on spines was verified. A CACTA transposon was found to be inserted into the first exon of the ws1 allele, knocking out its function and leading to the lack of spines in cultivated lettuce. All lettuce cultivars investigated have the nonfunctional ws1 gene, and a selection sweep was found at the WS1 locus, suggesting its important role in lettuce domestication. The expression levels of WS1 were associated with the density of spines among different accessions of wild lettuce. At least two independent loss-of-function mutations in the ws1 gene caused the loss of spines in wild lettuce. These findings provide new insights into the development of spines and facilitate the exploitation of wild genetic resources in future lettuce breeding programs.

Comprehensive identification and analysis of DUF640 genes associated with rice growth.

Protein domains with conserved amino acid sequences and uncharacterized functions are called domains of unknown function (DUF). The DUF640 gene family plays a crucial role in plant growth, particularly in light regulation, floral organ development, and fruit development. However, there exists a lack of systematic understanding of the evolutionary relationships and functional differentiation of DUF640 within the Oryza genus. In this study, 61 DUF640 genes were identified in the Oryza genus. The expression of DUF640s is induced by multiple hormonal stressors including abscisic acid (ABA), cytokinin (CK), ethylene (ETH), and indole-3-acetic acid (IAA). Specifically, OiDUF640-10 expression significantly increased after ETH treatment. Transgenic experiments showed that overexpressing OiDUF640-10 lines were sensitive to ETH, and seedling length was obstructed. Evolutionary analysis revealed differentiation of the OiDUF640-10 gene in O. sativa ssp. indica and japonica varieties, likely driven by natural selection during the domestication of cultivated rice. These results indicate that OiDUF640-10 plays a vital role in the regulation of rice seedling length.

miR396b/GRF6 module contributes to salt tolerance in rice.

Salinity, as one of the most challenging environmental factors restraining crop growth and yield, poses a severe threat to global food security. To address the rising food demand, it is urgent to develop crop varieties with enhanced yield and greater salt tolerance by delving into genes associated with salt tolerance and high-yield traits. MiR396b/GRF6 module has previously been demonstrated to increase rice yield by shaping the inflorescence architecture. In this study, we revealed that miR396b/GRF6 module can significantly improve salt tolerance of rice. In comparison with the wild type, the survival rate of MIM396 and OE-GRF6 transgenic lines increased by 48.0% and 74.4%, respectively. Concurrent with the increased salt tolerance, the transgenic plants exhibited reduced H2 O2 accumulation and elevated activities of ROS-scavenging enzymes (CAT, SOD and POD). Furthermore, we identified ZNF9, a negative regulator of rice salt tolerance, as directly binding to the promoter of miR396b to modulate the expression of miR396b/GRF6. Combined transcriptome and ChIP-seq analysis showed that MYB3R serves as the downstream target of miR396b/GRF6 in response to salt tolerance, and overexpression of MYB3R significantly enhanced salt tolerance. In conclusion, this study elucidated the potential mechanism underlying the response of the miR396b/GRF6 network to salt stress in rice. These findings offer a valuable genetic resource for the molecular breeding of high-yield rice varieties endowed with stronger salt tolerance.

Identification and characterization of BES1 genes involved in grain size development of Oryza sativa L.

BES1 (BRI1-EMS-SUPPRESSOR1) defines a unique class of plant-specific transcription factors that plays an essential role in response to Brassinosteroids (BRs) signal induction pathways. In our study, we conducted genome-wide scanning and comprehensive characterization of the BES1 gene family in rice and other eukaryotes, leading to valuable findings. Molecular docking experiments showed that all OsBES1 genes in rice could directly bind to BR small molecules. Among the identified genes, OsBES1-4 exhibited a remarkable response as it consistently showed induction upon exposure to various phytohormones after treatment. Further functional verification of OsBES1-4 revealed its impact on grain size. Overexpression of OsBES1-4 resulted in increased grain size, as confirmed by cytological observations showing an increase in cell length and cell number. Moreover, we identified that OsBES1-4 plays a role in rice grain size development by binding to the BR response element in the promoter region of the OsBZR1 gene. Evolutionary analysis indicated differentiation of OsBES1-4 between indica and japonica rice varieties, suggesting natural selection during the domestication process of cultivated rice. Therefore, we conclude that OsBES1-4 plays a crucial role in regulating rice grain size and has the potential to be an important target in rice breeding programs, and haplotype analysis found that all OsBES1 genes were associated with grain size development, either thousand-grain weight, grain length, or grain width. Overall, these findings suggest that the BES1 genes are involved in the regulation of grain size development in rice, and the utilization of SNPs in the OsBES1-4 gene promoter could be a favorable option for distinguishing indica and japonica.

A transposon-derived gene family regulates heading date in rice.

As a consequence of transposon domestication, transposon-derived proteins often acquire important biological functions. However, there have been limited studies on transposon-derived proteins in rice, and a systematic analysis of transposon-derived genes is lacking. Here, for the first time, we conducted a comprehensive analysis of the DDE_Tnp_4 (DDE) gene family, which originated from transposons but lost their transpositional ability and acquired new gene functions in Oryza species. A total of 58 DDE family genes, categorized into six groups, were identified in Oryza species, including 13 OsDDE genes in Oryza sativa ssp. japonica. Our analysis indicates that gene duplication events were not the primary mechanism behind the expansion of OsDDE genes in rice. Promoter cis-element analysis combined with haplotype analysis confirmed that OsDDEs regulate the heading date in rice. Specifically, OsDDE9 is a nuclear-localized protein expressed ubiquitously, which promotes heading date by regulating the expression of Ghd7 and Ehd1 under both short-day and long-day conditions. Single-nucleotide polymorphism (SNP) variations in the OsDDE9 promoter leads to changes in promoter activity, resulting in variations in heading dates. This study provides valuable insights into the molecular function and mechanism of the OsDDE genes.

Corrigendum: Comprehensive evolutionary analysis of growth-regulating factor gene family revealing the potential molecular basis under multiple hormonal stress in Gramineae crops.

[This corrects the article DOI: 10.3389/fpls.2023.1174955.].

Comprehensive evolutionary analysis of growth-regulating factor gene family revealing the potential molecular basis under multiple hormonal stress in Gramineae crops.

Growth-regulating factors (GRFs) are plant-specific transcription factors that contain two highly conserved QLQ and WRC domains, which control a range of biological functions, including leaf growth, floral organ development, and phytohormone signaling. However, knowledge of the evolutionary patterns and driving forces of GRFs in Gramineae crops is limited and poorly characterized. In this study, a total of 96 GRFs were identified from eight crops of Brachypodium distachyon, Hordeum vulgare, Oryza sativa L. ssp. indica, Oryza rufipogon, Oryza sativa L. ssp. japonica, Setaria italic, Sorghum bicolor and Zea mays. Based on their protein sequences, the GRFs were classified into three groups. Evolutionary analysis indicated that the whole-genome or segmental duplication plays an essential role in the GRFs expansion, and the GRFs were negatively selected during the evolution of Gramineae crops. The GRFs protein function as transcriptional activators with distinctive structural motifs in different groups. In addition, the expression of GRFs was induced under multiple hormonal stress, including IAA, BR, GA3, 6BA, ABA, and MeJ treatments. Specifically, OjGRF11 was significantly induced by IAA at 6 h after phytohormone treatment. Transgenic experiments showed that roots overexpressing OjGRF11 were more sensitive to IAA and affect root elongation. This study will broaden our insights into the origin and evolution of the GRF family in Gramineae crops and will facilitate further research on GRF function.

Creation of Elite Rice with High-Yield, Superior-Quality and High Resistance to Brown Planthopper Based on Molecular Design.

Breeding rice (Oryza sativa L.) with high yield, superior quality, desired grain shape and high resistance is the goal of breeding to meet the needs of current consumers. It is usually hard to combine multiple complex traits based on traditional breeding methods because they are frequently antagonistic to each other. However, molecular design breeding, as a novel breeding method, is an optional alternative to this challenge. To demonstrate molecular design breeding, 15 favorable genes from five parent lines were pyramided together to develop elite rice with high-yield, superior-quality, desired grain shape and high resistance to brown planthopper (BPH). The parental lines were 9311, the recurrent parent, carrying APO1, Ghd7, Ghd8 and Gn1a for high yield, GS3 and qSW5 for grain shape, and Wx and ALK for eating and cooking quality; 1880 with Gn8.1 for large panicles; Luo-Yu-Xiang carrying GW7 for grain shape and SBE3, SSIV2 and SSIII for eating and cooking quality; Luoyang6 with Bph6 and Luoyang9 with Bph9 for BPH resistance. After careful screening for the 15 targeted genes, desired phenotype and maximum genetic background from 9311, three molecular design lines with desired phenotypes, named as MD1 (Molecular design 1), MD2 and MD3 were developed. MD3 carried all 15 targeted genes, and MD1 and MD2 had 14 of the 15 targeted genes. Only SBE3 was not introgressed into MD1 and MD2 but this had minimal impact on the gel consistency and alkali spreading value. These newly bred lines exhibited higher yield potential, better grain quality with slender grains, low amylose content, high gel consistency and alkali spreading value, and higher BPH resistance compared to the parent 9311. In this study, we successfully created three novel rice lines with high yield, superior quality and improved BPH resistance by rational molecular design. Our results demonstrate molecular design is a powerful strategy to improve multiple complex traits and will provide a reference for the future commercial rice improvement.

Genome-Wide Identification and Analysis of the Metallothionein Genes in Oryza Genus.

Metallothionein (MT) proteins are low molecular mass, cysteine-rich, and metal-binding proteins that play an important role in maintaining metal homeostasis and stress response. However, the evolutionary relationships and functional differentiation of MT in the Oryza genus remain unclear. Here we identified 53 MT genes from six Oryza genera, including O. sativa ssp. japonica, O. rufipogon, O. sativa ssp. indica, O. nivara, O. glumaepatula, and O. barthii. The MT genes were clustered into four groups based on phylogenetic analysis. MT genes are unevenly distributed on chromosomes; almost half of the MT genes were clustered on chromosome 12, which may result from a fragment duplication containing the MT genes on chromosome 12. Five pairs of segmental duplication events and ten pairs of tandem duplication events were found in the rice MT family. The Ka/Ks values of the fifteen duplicated MT genes indicated that the duplicated MT genes were under a strong negative selection during evolution. Next, combining the promoter activity assay with gene expression analysis revealed different expression patterns of MT genes. In addition, the expression of OsMT genes was induced under different stresses, including NaCl, CdCl2, ABA, and MeJ treatments. Additionally, we found that OsMT genes were mainly located in chloroplasts. These results imply that OsMT genes play different roles in response to these stresses. All results provide important insights into the evolution of the MT gene family in the Oryza genus, and will be helpful to further study the function of MT genes.

Integrative Transcriptomic and Proteomic Analysis Reveals an Alternative Molecular Network of Glutamine Synthetase 2 Corresponding to Nitrogen Deficiency in Rice (Oryza sativa L.).

Nitrogen (N) is an essential nutrient for plant growth and development. The root system architecture is a highly regulated morphological system, which is sensitive to the availability of nutrients, such as N. Phenotypic characterization of roots from LY9348 (a rice variety with high nitrogen use efficiency (NUE)) treated with 0.725 mM NH4NO3 (1/4N) was remarkable, especially primary root (PR) elongation, which was the highest. A comprehensive analysis was performed for transcriptome and proteome profiling of LY9348 roots between 1/4N and 2.9 mM NH4NO3 (1N) treatments. The results indicated 3908 differential expression genes (DEGs; 2569 upregulated and 1339 downregulated) and 411 differential abundance proteins (DAPs; 192 upregulated and 219 downregulated). Among all DAPs in the proteome, glutamine synthetase (GS2), a chloroplastic ammonium assimilation protein, was the most upregulated protein identified. The unexpected concentration of GS2 from the shoot to the root in the 1/4N treatment indicated that the presence of an alternative pathway of N assimilation regulated by GS2 in LY9348 corresponded to the low N signal, which was supported by GS enzyme activity and glutamine/glutamate (Gln/Glu) contents analysis. In addition, N transporters (NRT2.1, NRT2.2, NRT2.3, NRT2.4, NAR2.1, AMT1.3, AMT1.2, and putative AMT3.3) and N assimilators (NR2, GS1;1, GS1;2, GS1;3, NADH-GOGAT2, and AS2) were significantly induced during the long-term N-deficiency response at the transcription level (14 days). Moreover, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that phenylpropanoid biosynthesis and glutathione metabolism were significantly modulated by N deficiency. Notably, many transcription factors and plant hormones were found to participate in root morphological adaptation. In conclusion, our study provides valuable information to further understand the response of rice roots to N-deficiency stress.

Quantitative Trait Loci Mapping of Mineral Element Contents in Brown Rice Using Backcross Inbred Lines Derived From Oryza longistaminata.

Mineral elements play an extremely important role in human health, and are worthy of study in rice grain. Wild rice is an important gene pool for rice improvement including grain yield, disease, and pest resistance as well as mineral elements. In this study, we identified 33 quantitative trait loci (QTL) for Fe, Zn, Se, Cd, Hg, and As contents in wild rice Oryza longistaminata. Of which, 29 QTLs were the first report, and 12 QTLs were overlapped to form five clusters as qSe1/qCd1 on chromosome 1, qCd4.2/qHg4 on chromosome 4, qFe5.2/qZn5.2 on chromosome 5, qFe9/qHg9.2/qAs9.2 on chromosome 9, and qCd10/qHg10 on chromosome 10. Importantly, qSe1/qCd1, can significantly improve the Se content while reduce the Cd content, and qFe5.2/qZn5.2 can significantly improve both the Fe and Zn contents, they were delimited to an interval about 53.8 Kb and 26.2 Kb, respectively. These QTLs detected from Oryza longistaminata not only establish the basis for subsequent gene cloning to decipher the genetic mechanism of mineral element accumulation, but also provide new genetic resource for rice quality improvement.

Quantitative Trait Locus Mapping of the Combining Ability for Yield-Related Traits in Wild Rice Oryza longistaminata.

In decades of hybrid rice breeding, the combining ability has been successfully used to evaluate excellent parental lines and predict heterosis. However, previous studies for the combining ability mainly focused on cultivated rice and rarely involved wild rice. In this study, for the first time, we identified 20 new quantitative trait loci (QTLs) for the combining ability in wild rice using a North Carolina II mating design. Among them, qGCA1, one of the major QTLs that can significantly improve the general combining ability of the plant height, spikelet number, and yield per plant, was delimited to an interval of about 72 kb on chromosome 1. qSCA8, another major QTL, which can significantly improve the specific combining ability of the seed-setting rate and yield per plant, was located in an interval of about 90 kb on chromosome 8. These QTLs discovered from wild rice will provide new ideas to explain the genetic mechanism of the combining ability and establish the basis for breeding of high-combining-ability rice.